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Abstract 

For single-antenna Gaussian interference channels, we re-formulate the problem of determining the Generalized 
Degrees of Freedom (GDoF) region achievable by treating interference as Gaussian noise (TIN) derived in [3] from a 
combinatorial optimization perspective. We show that the TIN power control problem can be cast into an assignment 
problem, such that the globally optimal power allocation variables can be obtained by well-known polynomial time 
algorithms (e.g., centralized Hungarian method or distributed Auction algorithm). Furthermore, the expression of 
the TIN-Achievable GDoF region (TINA region) can be substantially simplified with the aid of maximum weighted 
matchings. We also provide conditions under which the TINA region is a convex polytope that relax those in [3]. 
For these new conditions, together with a channel connectivity (i.e., interference topology) condition, we show TIN 
optimality for a new class of interference networks that is not included, nor includes, the class found in [3]. 

Building on the above insights, we consider the problem of joint link scheduling and power control in wireless 
networks, which has been widely studied as a basic physical layer mechanism for device-to-device (D2D) communi¬ 
cations. Inspired by the relaxed TIN channel strength condition as well as the assignment-based power allocation, 
we propose a low-complexity GDoF-based distributed link scheduling and power control mechanism (ITLinQ-b) 
that improves upon the ITLinQ scheme proposed in [4] and further improves over the heuristic approach known as 
FlashLinQ. It is demonstrated by simulation that ITLinQ-i- provides significant average network throughput gains over 
both ITLinQ and FlashLinQ, and yet still maintains the same level of implementation complexity. More notably, the 
energy efficiency of the newly proposed ITLinQ-i- is substantially larger than that of ITLinQ and FlashLinQ, which is 
desirable for D2D networks formed by battery-powered devices. 
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1. Introduction 

Power control and treating interference as Gaussian noise (TIN) is one of the most well-known, vastly employed, 
and yet most attractive interference management techniques, due to its low complexity, robustness to channel 
uncertainty, and to the fact that codes for the single-user Gaussian channel are well understood and efficiently 
implemented. Interestingly, it has also been shown that in some cases TIN is optimal or approximately optimal. For 
example, we know that TIN achieves the sum-capacity in the noisy regime of the two-user Gaussian interference 
channel [5]-[7]. In the general Ff-user single-antenna Gaussian interference channel, Geng et. al [3] have shown 
that, subject to a certain set of conditions on the channel strengths, TIN achieves the optimal Generalized Degrees 
of Freedom (GDoF) region, and achieves the capacity region to within a constant gap, independent of the channel 
coefficients and the signal-to-noise ratio (SNR). The TIN optimality condition found in [3] is simply expressed in 
words as the fact that, for each user (i.e., intended transmitter-receiver pair) the desired signal strength level is no less 
than the sum of maximum strengths of all interfering signals from the transmitter to the other (unintended) receivers, 
and to the receiver from the other (unintended) transmitters, when all signal strengths are expressed in log-scale 
(e.g., in dB). For future reference, we indicate this condition as the “GNAJ” condition, from the initials of the 
authors of [3]. Under the GNAJ condition, the TIN-Achievable GDoF region (briefly referred to as “TINA region”) 
is a convex polytope defined by the individual GDoF constraints and by the sum GDoF inequalities corresponding to 
all possible ordered subsets of users. With the aid of a combinatorial tool named potential graphs, the Ff-user TINA 
region was characterized in [3] by Em =2 (m)(^ -!)!?« (K - 1)! constraints. More recently, it has been also shown 
by Sun and Jafar in [8] that, by a series of transformations of linear programs, the sum-GDoF characterization can 
be translated into a minimum weighted matching problem in combinatorial optimization. As such, the sum-GDoF 
under the GNAJ condition can be characterized as disjoint cycles partition of the interference network. 

Such remarkable findings have inspired various related works, such as the TIN optimality of general X-channels 
[9], parallel interference networks [8], and compound interference networks [10]. In general, the TIN problem 
consists of two subproblems. Beyond the TINA region characterization, it is also important to find efficient methods 
to solve the TIN power control problem, that is, finding the (minimum) transmit powers that achieve a certain 
desired GDoF-tuple in the TINA region. The TIN power control problem has been open for a long time until a 
recent progress reported by Geng and Jafar in [10], where a simple yet elegant polynomial-time centralized iterative 
algorithm to find the globally optimal power allocation variables is provided. This centralized algorithm relies still 
on the representation by potential graphs. 

One may wonder if the potential graph representation is the only path to both TINA region characterization and 
TIN power control problems. Further, due to the distributed nature of interference channels, decentralized power 
allocation algorithms are more interesting, desirable and yet challenging. In addition, it is worth noting that the 
GNAJ condition was only proven to be sufficient. An interesting counter-example in [3] showed that there exist 
partially-connected (in the sense of channel strength levels) interference channels, such that TIN achieves the optimal 
GDoF region and yet the GNAJ condition is not satisfied. A natural question then arises as to whether there exists a 
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larger class of networks, including partially-connected ones, such that TIN is GDoF-optimal (i.e., TIN with power 
control still achieves the optimal GDoF region of the channel). These questions motivate this work. 

In this paper, the optimality of TINA is revisited. The TIN optimality problem was formulated in [3] by first 
eliminating power allocation variables using the potential theorem [11], to establish the TINA region in terms of 
GDoF variables only, and then by finding the optimal power allocation variables for a given GDoF-tuple in the 
TINA region [10]. In contrast, we re-formulate this problem in a reversed way, from a combinatorial optimization 
perspective [11]. Interestingly, by first casting power allocation into an assignment problem, the globally optimal 
power allocation variables corresponding to any feasible GDoF tuple in the TINA region can be found by solving 
the equivalent assignment problem in polynomial time, either in a centralized manner (e.g., Hungarian method 
[12], [13]) or in a distributed one (e.g.. Auction algorithm [14]). Inspired by the duality between the assignment 
and the maximum weighted matching problems in combinatorial optimization [15], we can express the TINA 
characterization in terms of a maximum weighted matching problem. In doing so, the TINA region is significantly 
simplified, requiring only 2^-1 constraints instead of (iT - 1)!. Interestingly, such a representation also offers 
an interpretation of the disjoint cycle partition in [8]. By this new formulation, we show that the TINA region 
is a convex polytope under a novel channel strength condition that relaxes the GNAJ condition in [3]. This new 
condition requires that the desired signal strength of each user is no less than the maximum difference between the 
sum strength of any pair of incoming/outgoing interference signals and the strength of the link between such a pair 
(all in dB scale). Furthermore, together with a connectivity condition, we are able to establish the optimality of 
TINA for a new class of networks. Such conditions are not included nor include the GNAJ condition [3]. 

Whereas fascinating from a conceptual point of view, how to translate these results into practice is also of great 
interest and practical importance to system designers. Device-to-Device (D2D) communication is expected to play 
an important role in future wireless communication systems (e.g., 5G), including applications such as car-to-car, 
machine-to-machine, proximity-based services, and multi-hop infrastructureless mesh networks. The physical layer 
of D2D systems is usually modeled as a Gaussian interference channel. Under the practical constraint of treating 
interference as Gaussian noise for the sake of complexity and robustness,^ a long-standing problem consists of 
controlling the power of the D2D links (transmit-receive pairs) in order to maximize the overall network throughput.^ 
The usual approach of guaranteeing a target signal-to-interference plus noise ratio (SINK) to each link turns out to 
yield an operating point that can be arbitrarily far from optimal. This is because some bottleneck links may impose 
too stringent constraints to the overall network. In contrast, much better network throughput can be achieved by 
selecting a subset of active links in each slot and allocating positive power only to these selected links [4], [18], 
[19]. By scheduling the subsets of active links over time, it is possible to achieve individual throughputs such that 
some network utility function is maximized. In turn, the shape of the network utility function determines the desired 

^From [16] we know that this condition is essentially equivalent to imposing the use of minimum distance decoding at each receiver. 

^Consistently with [17], we use the term “throughput” to indicate the time-averaged rate over a long sequence of scheduling time slots. In 
contrast, the instantaneous rate is the rate achieved in a single slot, for a given set of active users, i.e., links with positive transmit power. 
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fairness criterion (e.g., see [20], [21]). Link selection and scheduling has become the subject of intensive research. 
This problem is closely related to power control, since link selection corresponds to allocating zero or positive 
power to the transmitters. For a general D2D network, this problem is non-convex and, as a mater of fact, has a 
combinatorial nature. For example, a well-known power control method consists of replacing log(l + SINK) with 
log (SINK) in the user rate expression, and using Geometric Programming (GP) [22]. However, by neglecting the 
“1+” inside the “log” one has implicitly forced all links to use positive power, since assigning zero power to some 
links would drive the GP objective function to -oo. Instead, it is known that generally much better solutions can be 
found by first selecting a “good” subset of active links, and then allocating (positive) power only to the selected 
links. 

Various schemes for link selection have been proposed in the literature, e.g., [4], [18], [23]-[25] to name a few. 
For example, a large number of works is based on constructing an interference conflict graph [26], and then selecting 
maximal independent sets. These “maximal independent set scheduling” schemes are flawed by a fundamentally 
arbitrary choice of the threshold according to which two links are considered to be in conflict. Recently, a distributed 
link scheduling mechanism called FlashLinQ was proposed in [18], using a more dynamic link selection policy. 
Compared to those “maximal independent set scheduling” schemes, FlashLinQ dynamically takes both signal and 
interference strength into account. In FlashLinQ, links are ranked in priority order and considered one by one. 
Whether or not a link will be scheduled depends on whether this link does not cause/receive too much interference 
to/from links of higher priority that have already been selected (i.e., declared active). More recently, inspired by 
the GNAJ condition in [3], the authors in [4] proposed a new distributed link scheduling mechanism (referred to 
as “ITLinQ”) that provides significant sum throughput gains over FlashLinQ and yet maintains the same level 
of low-complexity. Instead of comparing the ratio of signal to interference strength of the new link with a fixed 
threshold as in FlashLinQ, ITLinQ compares the interference level caused to/received from existing links with an 
appropriately chosen exponent of the signal strength of the new link. It was verified by simulation in [4] that ITLinQ 
outperforms FlashLinQ with 28%-110% gains for a scenario where up to 4096 links can be scheduled. 

As a matter of fact, for general channel strength coefficients, the maximal subset of links satisfying the GNAJ 
condition may not lead to the maximal (weighted) sum throughput or sum GDoF. As will be demonstrated later, our 
relaxed channel strength conditions provide a larger convex polytope TINA region. This provides a generally larger 
subset of links on which power control can be applied, resulting in higher weighted sum GDoF. As a consequence, 
we are able to design a new distributed link scheduling and power control mechanism (named “ITLinQ+”), further 
fine-tuning the decision criterion of link selection. It is demonstrated by simulation that, without power control, 
ITLinQ-i- gains 5%-20% average sum throughput improvement over ITLinQ with 1024 links, at the expense of 
limited signaling overhead. When we also include power control, the average sum throughput is further enhanced, 
and more notably, the energy efficiency of ITLinQ+ is substantially improved (e.g., 50-100 times improvement for a 
10-link D2D network). In short, ITLinQ-i- improves the sum throughput performance and yet requires much less 
energy consumption, which is desirable for battery-powered D2D communications. Notice that achieving better or 
equal throughput with less energy consumption is not a contradiction here, since the network is operated in an 
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interference limited regime, such that rate is not immediately and obviously related to transmit power. 

This paper is organized as follows. In the next section, we present the system model of the general if-user 
Gaussian interference channels, followed by a summary of the main existing results of the approximate optimality 
of treating interference as Gaussian noise. In Section III, we reformulate the TIN problem from a combinatorial 
optimization perspective and we obtain a simplified description of the TINA region. By the simplified TINA region, 
we are able to identify a relaxed channel strength condition under which the general TINA region is a convex 
polytope. In Section IV, we consider the GDoF-based link scheduling and power control problem, offering a 
framework in this regard. Driven by this framework, the new decentralized link scheduling and power control 
mechanism named ITLinQ+ is proposed in Section V with detailed implementations. Section VI provides numerical 
results and comparisons with ITLinQ and FlashLinQ for some scenarios of D2D networks. We conclude the paper 
in Section VII. 

Notation : Throughout this paper, we define /C = {l,2,...,if}. Let A, A, and A represent a variable, a set, and 
a matrix, respectively. In addition, is the complementary set of A, and \A\ is the cardinality of the set A. Aij 
presents the ij-th entry of the matrix A, and A^ is the i-th row of A. As - ^5}, and As - Define 

A\a = {x\x G A,X ^ a} and Ai\A 2 - {x\x g Ai,x ^ A 2 }. Logarithms are in base 2. With a bit abuse of notation, 
k t i t j means k t i, i ^ j and k + j. 


11. System Model 


A. Channel Model 

We consider a if-user interference channel where both transmitters (Tx) and receivers (Rx) are equipped with a 
single antenna each. We shall refer to the j-th Tx-Rx pair as the j-th user. At Rx-j (VjfG/C = {l,...,if}), the 
received signal at the discrete-time instant t is given by 


K 


( 1 ) 


i=l 

where Xi{t) is the transmitted signal from Tx-i with power constraint E (|V^(f)p) < P^, hij is the channel coefficient 
between Tx-i and Rx-j, Zj{t) ~ CAi(0,1) is the (normalized) additive white Gaussian noise at Rx-j. Following [3], 
we translate the signal model in (1) into an equivalent GDoF-friendly form, given by 

Yj{t) = ^ + Z^{t) 


( 2 ) 


where Xi{t) = is the normalized transmitted signal with power constraint E(|Xj(t)p) < 1, \/and Oij 
are magnitude and phase of the channel coefficient between Tx-i and Rx-j, respectively, and the exponent aij is 
defined as the corresponding channel strength level 


log(max{l,|/ijjpPj}) 

logP 


( 3 ) 


where P > 1. Given a transmit power P^* at Tx-i with < 0, the Signal to Interference plus Noise Ratio (SINR) 


achieved by TIN at Rx-j is given by 


pOCjj+Vj 


pOCip 


We assume that the transmitters know channel strength levels 


perfectly for power control, and the receivers have access to both the magnitude and phase of channel coefficients. 
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B. Treating Interference as Noise 


We follow standard definitions for encoding/decoding functions and achievable rates. The individual achievable 
GDoF of message Wk is defined as dk - limp^oo where Rk is the achievable rate of user k. The (optimal) 
GDoF region V is the collection of all achievable GDoF-tuples (di, ^ 2 ,..., The TIN-Achievable GDoF 
(TINA) region defined in [3] is the set of all AT-tuples ((ii,d 2 ,... with components satisfying 


dj < max 


] 0, ajj + Vj - max{0, max( 


aij 



(4) 


for some assignment of the power allocation variables (ri, r 2 ,..., r^) ^ In the following, we denote the TINA 
by where the dependence on the specific network defined by {aij : i,j e JC} is clear from the context. 

From [3] we also know that the polyhedral TINA region is obtained by removing the positive part operator^ from 
the right-hand side of (4). Using the potential theorem [11], the authors of [3] are able to find a convex polytope 
form for the polyhedral TINA region for any subnetwork formed by a subset S ^ 1C and its associated desired and 
interfering links. We shall denote such poly tope by Since removing the positive part in the right-hand side 

of (4) restricts the GDoF region, then is achievable by switching off all users in = 1C\S and by using 

TIN for the users in S. We also denote by Vf the optimal GDoF region of the interference network, i.e., the region 
of GDoF-tuples achievable over any possible coding scheme (not restricted to TIN). 

The main results in [3] are summarized as below. 


Theorem 1. [ GNAT [3] ^ Consider a K-user single-antenna Gaussian interference channel with channel strengths 

{aij :ij eK}. 

1) For any subnetwork formed by users in S ^fC, can be described by^ 

0 < d/c < akk, V/c € 5, di = 0, Vi g 

m-l m-1 

k=0 k=0 

V ordered subsets ... ^ m g {2,..., |5|}. (5) 


2 ) The TINA region of the whole network is given by 


^TINA ^ pTINA 
5c/C 


( 6 ) 


3) If^kelC, 


akk > mdix{aik} + max{(A/cj}, (7) 

i:itk j-j^k 

then TIN is GDoF-optimal, i.e., Vf - T^tina _ ^tina whole region is a single convex polytope). 


^The positive part of x is max{0,a:}. 

^We use the term ordered subset to indicate that order matters, but elements are not repeated. For example, (1, 2, 3) and (1, 3, 2) are two 
rising such subsets for m - 3, but (1,2,2) is not valid, because it contains repeated elements. 
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Remark 1. It is easy to see that, for S = JC, there are in total Em =2 (m)(^ - 1)! {K - 1)! constraints in (5). 
Since the sum-GDoF dik does not depend on the order of the indices, for each unordered set of indices 

{io, • • • there are {m- 1)! inequalities, of which only one is relevant. However, finding which one is relevant 

involves, in general, extensive search, such that finding a general more compact form that eliminates redundant 
inequalities is non-trivial. 

III. TIN Problem Reformulation from a Combinatorial Perspective 

The expression of the TINA region in (5) involves a huge number of constraints, some of which are redundant. 
However, it is unclear which one is necessary and which one is required, such that we do not really know how 
to analyze this region. To make progress in this regard, we re-formulate the TIN problem of [3], [10] from a 
combinatorial optimization perspective. By casting the power allocation into an assignment problem, we find an 
alternative form for the TINA region via its dual - the maximum weighted matching problem [15]. Some basic 
definitions of weighted matching are recalled in Appendix VII-A. 

A. Casting Power Allocation into Assignment Problems 

In what follows, we consider a feasible GDoF tuple in for any user set 5 c /C, where ^ 

dj = ajj + Tj - max{0, max((T^j + r^)}, j g 5 (8) 

i:Hj 

given power allocation parameters {vj^j g 5}. By introducing two sets of auxiliary variables, namely, left labels 
{Vu ,} and right labels {y ^.} 



II 

1 

(9) 


t/„. = max{0, max(aij + r*)}, 

(10) 

the individual achievable GDoF can be rewritten as 



dj = (^jj-iyuj +yvj)- 

(11) 

Thus, for 5 c /C, the feasibility of a 

GDoF-tuple can be guaranteed by the minimization of the auxiliary variables 

sum: 

min Y,{yu,+yv^) 

{yujiVvj} 

(12a) 


s-t. Uuj + Vvj — f^jj ~ dj , Vj G S. 

(12b) 


In general, a given GDoF-tuple in pJina achieved by different assignments of the power control variables 

{vj : j e )C}. The componentwise minimum configuration corresponding to a given target GDoF tuple is referred to 

^Note that we consider di > 0, i e S. If di = 0, user pair i will be not activated and we simply remove it from without affecting others. 


as the globally optimal power control assignment. In this case, no users can reduce its transmit power while still 
achieving the same GDoF-tuple. Using the fact that, for all i t j, 

Vui + Vvj = -Ti + inax{0, max(Q;i/j + ri^)} (13) 

> -Vi + m.wK(aifj + Vi') (14) 

— ^ij •) ( 15 ) 

we have the following theorem that solves the GDoF-based power control problem for a given feasible GDoF-tuple. 

Theorem 2. For any (dj : j e )C) e a feasible power allocation assignment {vj^j e S) can be found by 

solving the following linear program: 

(AP) : min +yv^) (16a) 

{yuj,Vvj}j^S 

s.t. Vui + Vvj ^ aij, Vi + j (16b) 

Vuj + Vvj > ajj - dj,^j 6 S (16c) 

Vuj > 0, Vvj > 0, Vj e S (16d) 

where Vj = -puj, Vj g S. This linear program can be recognized as a dual formulation of an assignment problem 
[15], so that the unique glob ally-optimal power allocation can be found in polynomial time (e.g., 0{K^)) using 
e.g., the (centralized) Hungarian method [12], [13] or the (distributed) Auction algorithm [14]. 

Remark 2. The general solutions of (16) may not lead to the globally optimal power allocation. Due to the relation 
'^j = -Puj Puj ^ O 5 Pvj ^0 for all j, we conclude that the problem of finding the globally minimal power 
allocation is equivalent to that of finding the maximum left label equilibrium (cf, max Y.j Puj) or the minimum 

right label equilibrium (cf, mhiY.jPvj) to minimize the sum of the overall labels (i.e., mhiY.j{Puj ^Pvj)) (16). 

The Hungarian method has various implementations, but most of them are dedicated merely to the minimization of 
the overall sum of both left and right labels. As said, the maximal left label equilibrium to minimize the overall 
sum of labels is more relevant in our context. Fortunately, the Kuhn-Munkres algorithm [12], [13], a variant of the 
Hungarian method, offers such an equilibrium in solving the assignment problem. A detailed implementation of the 
Kuhn-Munkres algorithm with some parameters specified to fit our problem is relegated to Appendix VIl-B (see 
Algorithm 2). It is also worth noting that, according to the equality in (11), the optimal solution to (16) is achieved 
when the equality of (16c) holds. This observation is also added to Algorithm 2 as the new termination criterion. 
In particular, for an assignment problem with size K, the Kuhn-Munkres algorithm requires at most K rounds of 
iteration to converge to the optimal assignment solution. This fact can be also used to check the feasibility of GDoF 
tuples: if Algorithm 2 does not converge to the optimal solution for a given GDoF tuple within K iterations, then 
this GDoF tuple is infeasible. 

It is also worthwhile to mention that a distributed Auction algorithm, originally due to Demange, Gale, and 
Sotomayor [14], achieves the minimum right label equilibrium, whose values are element-wise smaller than any 
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Other feasible ones, leading to the global optimality of power allocation in a decentralized manner. A detailed 
implementation is presented in Section V-C (see Algorithm 1). 


B. TINA Region Representation 


In the following, starting from the power allocation solution of Theorem 2 and exploiting the duality between 
assignment and maximum weighted matching problems, we shall re-formulate the TINA region in a more useful 
and compact form. First, given a GDoF tuple ((ii,..., and channel strength level values {aij^i^j e K}, we 


define the following matrix associated with the assignment problem (16): 


- 


ai 


it j 


ajj dj^ i - j 


(17) 


By the duality theory in linear programming, we observe that the dual problem of (16) is given by 


max ^ Aijx{i,j), 

(18a) 

s.t. ^ x{i,j)<l, 

ieU:{i,j)eS 

(18b) 

Y, 

jeV:{i,j)eS 

(18c) 

x(i,j) e [0,1]. 

(18d) 


It is known that this linear program has integer-valued optimal solutions for bipartite graphs [27]. Since in our 
case the graph associated to the transmitters and receivers in S and corresponding intended and interfering links is 
bipartite by construction, then (18) coincides with a maximum weighted matching problem, obtained by replacing 
(18d) with x{i,j) e {0,1}. 

Next, due to the complementary slackness condition [15], an edge (i, j) belongs to the maximum-weight matching, 
i.e., = 1, if and only if + Pvj = A^j. For a given S ^ IC, due to (11) we have that a feasible GDoF-tuple 

implies equality in (16c). Hence, it follows that the set of feasible GDoF-tuples (i.e, the region coincides 

with the set of all (dj : j e S) for which the maximum matching solution of (18) is {(j, j),j e 5}. By the strong 
duality theorem, the minimum of sum of all left and right labels {puj^l/vj} in the primal problem (16) is equal to 
the maximum sum weights of all matchings in the dual problem (18). This is the key observation that enables us to 
provide a more compact form for the TINA region. 

We construct a weighted full-connected bipartite graph 0 = (/C,/C,/Cx/C), where the weight a'- is specified as 


a. 


aij, itj 

0, i=j 


(19) 


For any 5 ^ /C, we define the subgraph ^[5] = (5,5,5 x 5) with weights {a[j e 5}. By the observation above, 
the sum of {Ajj : j g 5} must be no less than w(A4) for any matching A4 of ^[5]. Hence, we can write 


Y^{oLjj - dj) > maxu;(Af5) = w{M%), 


( 20 ) 
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where Ai% is the matching of G[S] with the maximum weight. This yields the following result. 

Theorem 3. Consider a K-user single-antenna Gaussian interference channel with channel strengths {aij : i^j e 1C}. 
For any user subset S ^ 1C, pTiNA ^ 

dk>0, ykeS, di=0,yieS^ 1 

\ ( 21 ) 

T.k€S'dk ^ T.k€S'^kk - J 

where = 0 |5'| = 1. This simplified representation is equivalent to the expression in (5). 

Proof: See Appendix VII-C. ■ 

Remark 3. For individual users, i.e., |5'| = 1, we have individual GDoF constraints, i.e., dk < (Xkk- Using Theorem 

3 into (6), we find that we need only 2^-1 non-trivial inequalities, one for each non-trivial subset of 1C, to describe 
the K-user TINA region which is significantly less than {K - 1)\ in [3]. 

Remark 4. [8, Theorem 3] states that if the GNAT condition (7) is satisfied, the sum-GDoF is equal to the best 
cyclic partition bound. Explicitly, the best cyclic partition of a user subset S is a partition S - {5i,... ,5^} ^ 
satisfying 

w(M*s) = Y,w(M*sJ- ( 22 ) 

2=1 


= (4 : fc e /C) : 



1.5 

(a) 


1.5 

0.5 

0.1 

0.2 

0.4 

0.5 

1 

0.5 

0.8 


( b ) 



Fig. 1: (a) A 3-user interference channel, and (b) the input weight matrix of Hungarian method for the GDoF tuple (0.5, 0.6, 0.7). 


Example 1. We consider the example in [10, Fig. 8] to show the efficiency of our formulation, as shown in Fig. 1(a). 
According to Theorem 3, the TINA GDoF region is immediately given as 

^{™ 3 } = {{dud2,d3) ■■ 0 < 4 < 2,0 < 4 < 1,0 < 4 < 1 - 5 , 
d\ 6^2 ^ 2.3,6^2 + d^ < 1.5 
d\ + 6^3 ^ 2.4, d\ + d2 + d^ ^ 2.5j', 

^Recall that {<Si,..., <Sp} is a partition of <S if n <5^ = 0 Vz ^ j and S = <3i. 
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which is identical to the expression found in [10]. In order to solve the power allocation for a given GDoF-tuple 
(say (0.5,0.6,0.7) in this case), we take the weight matrix in Fig. 1(b) as the input of the Kuhn-Munkres algorithm 
(see Algorithm 2 in Appendix VII-B) and we obtain: 

2/«i = 1-2, Vu-, = 0.4,2/„3 = 0.7,= 0.3, = 0, = 0.1. 

Thus, the globally optimal power allocation assignment is ri = -1.2, r 2 = -0.4, rs = -0.7, which coincide with what 
found in [10]. The details are relegated to Appendix VII-B. 

Clearly, to start Algorithm 2, yuj and yy. are initialized respectively with the maximum value of the j-ih row of 
A and 0. Following the procedure in Algorithm 2, we gradually decrease yuj and increase yy. to make sure the 
constraints in (16) satisfied. Note that Vj = -yuj is increasing during this procedure. Once we find one solution, it 
will be the global optimum assignment, because it is impossible to decrease rj (correspondingly increase yu^) and 
find another solution in the region that we have already explored. 0 

C. A New TIN Optimality Condition 

Besides the reduction of the number of inequalities, this new formulation enables us to identify a relaxed channel 
strength condition such that the TINA region is a convex polytope. 

Theorem 4. Consider a K-user single-antenna Gaussian interference channel with channel strengths {aij : i,j e 1C}. 

If 

akk > max {aik + akj - J, V ke 1C, (23) 

id- id*k 

where ol\- is defined in (19), then is monotonically non-decreasing with respect to S, i.e., if Si ^ S 2 ^ 1C 

then c Also, t^tina _ pTiNA ^ convex poly tope. 

Proof: See Appendix VII-D. ■ 

Remark 5. The newly found channel strength condition is a relaxed version of the GNAT condition (7), because 
aC is non-negative such that if (7) is satisfied, then (23) is satisfied automatically. When i = j + k, (23) reduces 
to akk ^ OLik + aki, V k,i, s.t. k ^ i. When i t j t k, it reduces to akk + o^ij > aik + OLkj, for all k t i, k + j 
and i ^ j. As such, in a network where the condition (23) is satisfied for all users, for a tuple {i,j,k) with 
i t j t k, we have either akk ^ + otkj or aik + o^kj > o^kk ^ + o^kj - otij. We conclude that, for a user k, if 

o^ik ^kj ^ Oikk — ^ik o^kj ~ OLij, then akk — oiki + or,jk- It is because, if both aik o^kj ^ ^/c/c ond aki oijk > oikk 
are satisfied, it leads to 2akk < o^ik + OLki + OLkj + OLjk, which conflicts with the fact that akk ^ OLki + cxik for all i + k. 

In view of the fact that (23) is a relaxation of (7), the corresponding TINA region (although a convex polytope) 
is not generally optimal. This is because when (23) holds but (7) does not, the converse argument to prove GDoF- 
optimality does not go through. However, we can exhibit a class of networks different from the class identified in 
[3], for which the TINA region is GDoF-optimal. This is a special class of partially connected interference channels 
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satisfying a topological condition given below. Interestingly, this class of networks is not included nor includes the 
class defined by the GNAJ condition. The converse proof follows the approach in [3] and is presented in Appendix. 

Theorem 5. Consider a K-user single-antenna Gaussian interference channel with channel strengths {aij : i^j e K,}. 
Assume that (23) holds and, in addition, that for every S with |5| > 2, and the corresponding fully connected 
weighted subgraph ^[5] = (5,5,5 x 5) with weights {aC 5}, 

3 (i, j) £ M*s, s.t. aij = 0. (24) 

Then, TZ* = 

Proof: See Appendix VII-E. ■ 

Remark 6. As the maximum weighted matching may not he unique. Theorem 5 holds as long as (24) holds for 
any one of the maximum matchings. Condition (24) allows us to establish the optimality of TINA since, under this 
condition, we can prove that the converse is tight. This, however, is only a sufficient condition and there might be a 
larger class of networks, including both the subclass defined by Theorem 1 and the one defined by Theorem 5, for 
which TIN is GDoF-optimal. 

Example 2. We illustrate the relaxed channel strength condition by the example in Fig. 2. It is easy to verify that the 
condition (23) holds for the entire network, while the original GNAJ condition (7) does not hold for users 1 and 2. 
Note that M"" = {(1,3), (2,1), (3,2)} is a (non-unique) maximum weighted matching and contains = ais = 0, 
such that also condition (24) holds. Thus, from Theorems 4 and 3, the TINA region of this network is the poly tope 
defined by: 

^{7,2,3} = {{dud2,di ): 0 < di < 1, Vi 6 {1,2,3} 

di + d2< 1.1, d2+ d^< 1.3 

di + 6^3 < 1 . 2 , di + 6^2 + d^ < 1 . 8 }. 

0 

Remark 7. A subclass of network topologies for which (24) holds is the class of networks that have no perfect 
matchings in any unweighted subgraph of Q with zero-weight edges removed. A bipartite graph has no perfect 
matchings if Hall’s condition does not hold [27]. The so-called triangular networks in [28] belong to this category. 

IV. A GDoF-based Fink Scheduling and Power Control Framework 

In this section we capitalize on the insight about the TINA region obtained before, in order to develop a framework 
for link scheduling and power control in Gaussian AT-users interference channels with the constraint that receivers 
treat interference as (Gaussian) noise. As anticipated in Section I, this finds applications in practical interference 
management of D2D networks, where devices communicate directly to their intended destinations sharing the same 
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Fig. 2: (a) A 3-user IC where TINA is a convex polytope. The value associated with each link represents the channel strength 
level aij and the missing links correspond to aij = 0. (b) The links where the relaxed channel strength condition (23) is 
satisfied while the GNAJ condition does not hold (marked in blue and purple). 


channel bandwidth. In general, the goal is to activate simultaneously a subset of links (i.e., transmitter-receiver 
pairs) with nonzero transmit power, aiming at maximizing some desired system utility function. The classical power 
control problem (e.g., as formulated in [29]-[31]) finds the componentwise minimum transmit power vector that 
achieves given target SINRs at the receivers, when such target SINRs are feasible. However, this approach does not 
take into account that in modern TDM A systems the links may not be active in all scheduling slots. In contrast, 
by selecting a subset of links on each slot (scheduling), higher user throughput (i.e., time-averaged rate) can be 
achieved. As anticipated in Section I, a direct application of GP [22] also does not solve the scheduling problem, 
since implicitly all links must be allocated positive power. Intuitively, these approaches work well when SINRs 
significantly larger than 1 (0 dB) can be achieved for all the K links. 

A general scheduling framework is provided by considering the user throughputs = limt^oo | Er=i 
where Rkir) indicates the rate achieved by link k during scheduling slot r. Let ... ,Tx) denote a concave 

componentwise non-decreasing Network Utility Function of the user throughputs, and let T denote the achievable 
throughput region of the system. Then, a general Network Utility Maximization (NUM) problem is given as 


(NUM): max ..., T^) 

(25a) 

s.t. (Ti,... ,Tk) e T. 

(25b) 

In our case, we shah consider a GDoF criterion and replace with dk = limp^oo 

Through an immediate 

time-sharing argument, we have that the achievable region of throughput-GDoF is the convex hull of denoted 

by convT^^^^^. In general, T^tina union of convex polytopes (see Theorem 1), such that it is not generally 

convex. However, when (23) in Theorem 4 holds, then t^tina _ _ ^tina Using the GDoF criterion. 

the corresponding NUM problem becomes 


(NUM-GDoF): max U(di,...,dK) 

(26a) 

s.t. (di,..., dx) ^ convT^^^^^. 

(26b) 


It turns out that the above problem can be solved by iterating over time (i.e., over the scheduling slot) a sequence of 
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“instantaneous” subproblems. The following result is quite standard and follows as corollary of the general theory 
developed for example in [21], [32]-[34] (and references therein), and shall be stated without proof here for the 
sake of space limitation. 


Theorem 6. Consider a K-user single-antenna Gaussian interference channel with channel strengths {aij 

and corresponding TINA region For a sequence of scheduling slots indexed by t = 1,2,3,..., consider the 

following iterative procedure: 

1) Initialize weights Wk{l) = 1 for all k e 1C. 

2) For f = 1,2,..., repeat the following two steps: 

• Compute the GDoF-tuple ..., (t)) solution of the max weighted sum GDoF problem 


{SUM - GDoF) : max ^k{t)dk 

keK 

s.t. 


(27a) 

(27b) 


• Update the weights according to 

Wk{t + l) =max{0,u;/e(t) -4(t) + a^(t)}, 

where ... ,a^(t)) is the solution of the convex optimization problem 

max VU{ai,...,aK)-YjWk{t)ak 
/C€/C 

S.t. (ui, . . . , Ox ) ^ [O? -^max] i 


(28) 


(29a) 

(29b) 


where V and ^max 
Then, for sufficiently large A 


> 0 are control parameters of the algorithm. 
max have that 


lim f/ I - 

t-^OO 




>U(d,,...,di,)-^, 


(30) 


—^ ^ 

where {d^,... ^dj^) is the solution of the NUM-GDoF problem (26), and k, is a constant that depends on the system 
parameters but is independent of V. Hence, the above iterative scheduling algorithm can approach the optimal 
value of (26) by any desired accuracy. 


Remark 8. It is also possible to show that the time over which the limit in (30) is closely approached grows as 
0{V). Therefore, in practice there is a tradeoff between how close we can approach the optimal network utility 
function value, and how quickly the scheduling algorithm converges. Nevertheless, here we are not concerned with 
this problem, and we use Theorem 6 as a general tool to translate a NUM problem in terms of the long-term time 
averaged rates (or GDoF, in our case) into a sequence of ''instantaneous” (i.e., to be solved at each scheduling 
slot) max weighted sum GDoF problem. 

It follows that, from now on, we shall be concerned with solving the max weighted sum GDoF problem (27) 
for an arbitrary set of weights (rci,... ,wk)- The “power control” aspect of the problem reside in the fact that. 
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when a solution , (i^(t)) of (27) is found, we must also provide the powers at which the links have to 

transmit in order to realize such GDoF point in the TINA region. As anticipated in Section I, such transmit powers 
are generally not unique, and in this case we aim at finding the globally optimal power control assignment for the 
desired GDoF-tuple. 

In what follows, we first introduce the exact GDoF-based solution, and then subsequently simplify it until we 
could obtain an approximation with polynomial-time complexity. In turns, the exact or approximate solver of 
(27) can be plugged into the iterative scheduling algorithm of Theorem 6 in order to obtain a scheme that works 
for any suitable network utility function. For example, if throughput max-min fairness is desired, we can choose 
f/((ii,..., dx) = min/c dk. Instead, if proportional fairness is desired, we can choose f/((ii,..., d^) = E/c€/c log(<^/c)- 

A. Exact Joint Solution is Hard 

We re-write (27) more conveniently in the form: 


max 

X! "^kdk 

(31a) 

{dk} 

keK. 


S.t. 

5c/C 

(31b) 


which can be categorized as an instance of Disjunctive Programming (DP) [35]. The union involves 2^-1 nontrivial 
polyhedra, and is described by 2^^^ -1 linear inequalities. As mentioned earlier, the union is not necessarily 

leading to a convex polytope, and thus the problem is not a convex optimization problem in general. Nevertheless, 
it can be transformed to an equivalent convex optimization problem by replacing with its convex hull 

Q = conv (U^sic 

The full description of Q may require an exponential number of inequalities, yet Q has a compact representation 
in a higher-dimensional space. The so-called lift-and-project cutting plane method [36], [37] can be employed to 
offer an exact solution to this problem. The principle consists of three steps: (1) lift the subspace spanned by GDoF 
tuples into a higher-dimensional space by introducing some auxiliary variables, (2) obtain the compact representation 
in the form of a set of lift-and-project cutting planes, and (3) project the compact representation onto the original 
GDoF spanned subspace. These cutting planes are valid for the closure of the convex hull Q, and can be generated 
by solving cutting generating linear programs derived from the higher dimensional representation (see [36], [37] and 
references therein). 

Once we obtain the GDoF tuple maximizing the weighted sum-GDoF, the second step is to use either the potential 
graph based centralized power allocation algorithm found in [3] or the assignment problem inspired algorithms (e.g., 
centralized Hungarian method or distributed Auction algorithm) presented in this paper to find the globally optimal 
power allocation parameters. This is illustrated in Fig. 3. 

Unfortunately, the cutting generating linear programs still involve exponential number of constraints. This fact 
prohibits the application of this exact solution via disjunctive programming for network of practical size (e.g., a few 
tens to a few hundreds of D2D links). As such, reasonable approximation and heuristic approaches are desirable, 
although the global optimality is not guaranteed. 
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Fig. 3: The exact solution, where the disjunctive program identifies the optimal GDoF tuple d in which the links in the subset S 
with positive GDoF values are scheduled, and the assignment problem outputs the power allocation vector r. 


B. Separated Link Scheduling and Power Control 

In view of the complexity of the exact solution, we resort to separated link scheduling and power control. We 
can first select heuristically a subset of links S whose TINA region contains a GDoF tuple that leads to a 

reasonably large weighted sum-GDoF for the given user profiles (i.e., weights). Then, we can find this GDoF tuple 
d by optimizing the linear program (i.e., the weighted sum-GDoF) within the TINA region Finally, given 

this GDoF tuple d, we can use the assignment inspired algorithms to obtain the globally optimal power allocation 
parameters r. This procedure is illustrated in Fig. 4. 


JC 
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Problem 



Fig. 4: The separated solution with a heuristic link scheduling and assignment-inspired power control, where the optimal GDoF 
tuple is obtained via linear programming over the selected convex polytope. 


1) Link Scheduling: Recall that in the disjunctive programming (31), the objective function can be regarded as a 
moving hyperplane, and the constraint is the union of polyhedra (convex polytopes). A locally optimal weighted 
sum-GDoF solution is met when the hyperplane touches one of the vertices as a tangent plane to one convex 
polytope. The vertices of the largest convex polytope meet such hyperplanes with high probability, such that the 
weighted sum-GDoF can be maximized with high probability in such a polytope. As such, a heuristic link scheduling 
is to select the largest subset of links whose corresponding TINA region is the largest polyhedron (convex polytope) 
among all polyhedra. 

Similarly to the so-called information theoretic independent sets (ITIS) introduced in [4], we define an independent 
set according to channel strength levels of our new TIN condition in Theorem 4. 

Definition 1. A user subset S is called an improved information theoretic independent set (referred to hereafter as 
‘TTIS+”), if for any link k e S 

ttkk > max {aik + akj - a', }. (32) 

i,jsS\{k} 

Remark 9. The improvement is in the sense that as long as a subset of users forms an ITIS, it automatically forms 
an ITIS+. The reverse claim is not true. As such, an ITIS+ contains no less links than an ITIS, and thus the TINA 
region of the former includes that of the latter. 
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The TINA regions of the selected users in ITIS+ are convex poly topes. The largest convex poly tope corresponds 
to the links that form the largest ISIT+. As such, link scheduling turns out to be finding the largest ISIT+. Although 
the GDoF optimality does not always hold in the largest ISIT+, it offers the potential to achieve a larger sum 
throughput than that obtained by ITIS. 

2) Power Control: As said earlier, we can first find the optimal GDoF tuple by maximizing the weighted 
sum-GDoF for the given set of scheduled links S over the convex polytope 77 >tina such a GDoF 

tuple, we can apply assignment inspired algorithms to obtain the globally optimal power allocation parameters. 

For a given ITIS+ S and the associated user profiles (e.g., weights of individual rates), the GDoF tuple with the 
maximum weighted sum GDoF can be identified by solving the following linear problem: 


max 

YjWidi 

(33a) 

di 

i&S 


S.t. 


(33b) 


For the linear programs, we have polynomial-time complexity algorithms in terms of the numbers of variables and 
constraints, such as simplex method, interior-point method. As there is usually an exponential number of constraints 
(i.e., - 1) in thus the complexity solving this linear program is still exponential with respect to |5|. 

The problem as to how to reduce the size of constraints by exploiting the special structure of linear program is an 
interesting problem yet beyond the scope of this paper. 

C. Replacing Linear Programming by Geometric Programming 

The exponential-time complexity of the linear problem in (33) prohibits its application in large D2D networks. To 
bypass this challenge, we replace it by a geometric program, which can be solved with polynomial-time complexity, 
together with a recalculation of GDoF. 

As illustrated in Fig. 5, given a selected subset of links 5, Linear Programming (LP) is replaced by Geometric 
Programming [22] (GP) to maximize the weighted sum rate at high SNR, and we obtain the corresponding power 
allocation vector r'. Basically, this power allocation r' achieves the maximum weighted sum throughput yet leads 
by no mean to the minimum power allocation. To obtain the minimum power allocation, we proceed further and 
regenerate the GDoF tuple d from the power allocation r' according to (8). Finally, assignment-inspired algorithms 
are utilized to obtain the minimum (globally optimal) power allocation parameters r. Admittedly, GP can be solely 
used for power control. Yet, by concatenating with the assignment-inspired algorithms, it lends itself to a lower 
power consumption, i.e., r' > r. 

In what follows, we establish the equivalence of such a replacement. Given the user profiles {wi^i e S} of a. 
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Fig. 5: The separated solution with a heuristic link scheduling and assignment-inspired power control, where the optimal GDoF 

tuple is obtained instead via geometric programming. 


selected use subset S, we initially aim at solving the following optimization problem: 


max 

Y,Wi\og{i + smRi) 

(34a) 

{P} 

i£S 



GiiPi 


s.t. 

SINRi = ^ ^ „ 

(34b) 


1 + Y,j^i GjiPj 



0 < Pi <1 

(34c) 

where Gij = \hij\^ is the channel gain between Tx-i and Rx-j. Let us introduce an auxiliary variable ti 

1 

SlNRi 


where U is a posynomial function of {Pi^i e S}. Thus, the optimization problem at high SNR can be approximated 
to 

(GP): min (35a) 

ts 

1 + GjiPj 

s.t. - ^ ^ ^ <ti, Vie<S (35b) 

GiiPi 

0<Pi<l (35c) 

which is a geometric program with regard to {Pi^UG ^ 5}. A brief description of GP can be found in Appendix. 
The equivalence between the LP in (33) and the GP in (35) is due to the following proposition. 

Proposition 1. Given a user subset S and the profile {wifi e 5}, if the condition in Theorem 4 holds, the power 
allocation {Pifi e S} by GP in (35) is equivalent to LP in (33), in the sense that both approaches achieve the same 
weighted sum rate at high SNR. 

Proof: Let Gij = max{l, and P^ = E (|X^(t)p) = P^" < 1. Substituting them into GP, we have 


min 

Yitr 

(36a) 


i€S 


S.t. 

1 -1- V T^Gtji+rj 

<ti 

pa-ii+Ti 

(36b) 


- oo < < 0 


(36c) 
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By replacing ti with P we have an equivalent formulation: 


max 

y,widi 

(37a) 

{ri,di} 

i€S 


S.t. 

^ _l_ ^ ^ j p~di+0'ii+ri 

(37h) 


jti 



Vi <0, \fi e S 

(37c) 


Note that the log-sum-exp function can be rewritten as 


log(l + ^P^^^^^O = (maxzi)logP + log^P"^"^^"^'4 

jti ^ j 


where = 


logE.P 


logP 


and 


= (maxP + e^) logP 

n ^ 


(38) 


1 Oiji+rj, j + i 

I 0, j = i 


(39) 


It is easily verified that 0 < e < , and thus the second term in RHS is always hounded within [0, logi^]. Thus, 

we can rewrite the linear program as: 


max 

y,widi 

(40a) 

{ri,di} 

us 


S.t. 

di < an + ri- max{0, (aji + rj)}-ei 

jtl 

(40b) 


di > 0, < 0, Vi e 5 

(40c) 


At high SNR (P ^ oo), ^ 0. It is not hard to verify that the feasible region of {di : i e S) in the above linear 
program is exactly the one in (33) formulated by taking GDoF metric into account. This completes the proof. ■ 


Remark 10. As mentioned in [10], for any rate tuple, there only exists a unique locally optimal power vector, which 
is also globally optimal (i.e., element-wise minimal), while there are multiple locally optimal power vectors for a 
GDoF tuple. That being said, multiple locally optimal power vectors lead to the same GDoF tuple, but only one 
locally optimal power vector is globally optimal, which can be obtained by assignment-inspired algorithms or the 
algorithms in [10]. The concatenation of the geometric program and assignment-inspired algorithms offers maximal 
sum throughput at high SNR as well as minimal power consumption. Note that they are not contradict, since the 
network is operated in an interference limited regime such that the sum rate is not immediately and obviously 
related to transmit power. 


V. ITLinQ-f: a Decentralized Implementation in D2D Communications 

In D2D communications, smart devices are distributively located such that they have to make their decisions 
uncoordinatedly. As such, a decentralized implementation of the link scheduling and power control framework 
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in Section IV is desirable and of great interest to system designers. FlashLinQ [18] and ITLinQ [4] are two 
low-complexity distributed link scheduling mechanisms for D2D networks with reasonable signaling overhead. 

In this section, we propose a decentralized mechanism (referred to as “ITLinQ-i-”) for the link scheduling and 
power control framework aforementioned in Fig. 5. Our proposed ITLinQ+ mechanism consists of three ingredients: 
(1) a decentralized implementation of link scheduling to find the largest ITIS+, (2) a decentralized GP implementation 
to find the power allocation vector and the corresponding GDoF tuple, and (3) a distributed Auction algorithm to 
solve the assignment problem and yield the globally nfinimal power allocation parameters. These ingredients are 
detailed in the following subsections. 


A. Decentralized Implementation of Link Scheduling 


To figure out the largest ITIS+, it requires to coordinate all devices and enumerate all possible combinations, which 
is expensive to coordinate and difficult to enumerate. A decentralized link scheduling criterion with a greedy-flavor 
independent sets selection, a reasonable amount of signaling, and a comparable complexity to FlashLinQ and ITLinQ, 
is highly demanded. 

The decentralized link scheduling of ITLinQ+ is comprised of two phases: link scheduling and signaling. We first 
present the link scheduling and then figure out the information that needs to pass. 

Link Scheduling Phase: We first order the links according to their priorities (e.g., weights of throughput). To 
initialize, the links with highest priority will be first selected, and new links will be added to them later. At a certain 
point, suppose we have selected some links in 5 = {A, • • •, Whether or not a new link ik is suitable to be 

scheduled depends on the following conditions. For the convenience of comparison with FlashLinQ and ITLinQ, we 
adopt the notations SNR/. = and INR^^ = ^ where the noise power is normalized. 

• At Tx-ik, check if the following condition is satisfied: 


SNR? > 


INR,. 


y j < k, ij e S 


(41) 


where 7^,7 e [0,1] are design parameters, and mins</c,s^j{INR^^i^.} is the least channel strength level of 
incoming interfering links of Rx-ij . 

At Rx-ik, check if the following condition is satisfied: 


SNR? > 


INR, 


'ij 


y j <k, ij e S 


(42) 


where mins</.^s^j{INRi^. 4 ^} is the least channel strength level of outgoing interfering links of Tx-ij. 

If these two conditions are satisfied, then this new link ik can be scheduled, i.e., S ^ Su{ik}. Note that the minimum 
value of INR is initialized to be 1 such that the second link 72 is selected if SNR?^ > max{INR ^^^2 5 } i^ 

satisfied. 

Link Signaling Phase: The potential performance improvement of ITLinQ+ over FlashLinQ and ITLinQ is at the 
expense of additional signaling among devices. Before the link scheduling phase, we have two rounds of signaling 
to inform transmitters and receivers the channel strength information, as did in FlashLinQ [18] and ITLinQ [4]. 
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• In the first round, similarly to FlashLinQ, the transmitted signals from all transmitters are sent with full power 
P in different frequency bands, such that each receiver is able to estimate the channel strength of the channels 
to which it is connected. 

• In the second round, similarly to ITLinQ, pilot signals are sent from all receivers with full power P from each 
receiver in different frequency bands. As such, the transmitters estimate the channel strength of the channels to 
which it is connected. 

At the end of this procedure, for any link k, the local channel strength information {INR/.^, Vi} and SNR/^ are 
accessible at transmitter k, and {INRj/e, V j} and SNR/^ at receiver k. 

Another additional signaling cost happens at the end of each successful link selection. The transmitter and 
receiver ij (ij e S) have to inform the next being checked links the minimum interfering channel strength (i.e., 
min^^e 5 ^ 5 ^j{INR^^^^} and min^^e 5 ^s^j{INR^^^^.}, respectively). This signaling can be done similarly as the above 
signaling procedure, or by broadcasting. 

To reduce this additional signaling overhead, we can replace both mmi^^s,stj{^^^isij} 
by min^^g^ , INR^^.^^}. The next being checked link ij^ is required to compare the minimum value among 

all cross links in S with the links to which it is connected (i.e., INR^^^^., Vij e S). If the minimum value does not 
change with the newly selected link, then transmitter or receiver does not need to do anything. Otherwise, link 
ik has to inform the next being checked links this updated minimum value. Finally, the decision criterion of this 
decentralized link scheduling can be stated in words: A link is scheduled if the interference caused to/received from 
the already selected higher-priority links is smaller than the product of the signal strength with an exponent p of this 
link and the signal strength with an exponent 7 of the weakest interfering link among the already selected links. 

B. Decentralized GP Implementation to Find the Optimal GDoF-tuple 

In what follows, we will consider a distributed implementation of GP [22] where Rx-i has access to the local 
knowledgejr'• = aji -\- rj}j only, i.e., the exponent of received signal power from Tx-j. From (34), we formulate a 
linear program taking into account the local knowledge 


max 

y + rj - max{0, r'J) 

(43a) 

ri 

lis 


S.t. 

i’ji = aji+rj,yj eS\{i}, Vi 

(43b) 


Vj < 0, Vj 6 S 

(43c) 


This can be solved by alternating optimization. First, given we obtain the optimal solution of local 

variables and {r'•, j e 5\{i}}. Then, given the updated local variables and {r'-, j € 5\{i}}, we update 
again. Keep doing this until it converges. 

Introducing Lagrange multipliers only for the coupling constraint, we form a partial Lagrangian 

L = -Y,Wi {an +ri- max{0, rji})+'Z Y, i^'ji “ “ ’’i ) 

i€S nSj^S\{i} 


(44) 
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and thus, each user only has to take care of its local partial Lagrangian term, given by 

Liin, {r'ji}j^i;{jji}) = -Wiiau + ri- max{0,^ -yjir'ji - ^ -yij n . (45) 

The minimization of partial Lagrangian can be done locally by each user in parallel with respect to the primal 
local variables and {rL, j g 5\{i}} given the knowledge of {'yji}j. 

The dual variable {jji}j can be obtained by solving the dual problem 

max ({Tji}i) - E ; {ij*})• (46) 

A simple solution of {'yji}j is to update iteratively with the updating rule in t iteration being: 

lji(t + 1) = (47) 

where = aji-\-rj{t) is the estimation of received signal power exponent, and 6{t) is a carefully chosen stepsize. 

Signaling of {'yji}j is needed in each iteration, and the reduction of such signaling overhead can be similarly done 
as in [22]. 

C Decentralized Auction Algorithm for Power Allocation 

The Auction algorithm is an iterative procedure to determine the optimal assignment of a number of products to 
a number of potential buyers fulfilling their own best interests. It mimicks the sales auction in the business activities 
in which bids are compared in multiple rounds to make the best offer to the products, with each product going to the 
highest bidder. The Auction algorithm is an efficient way to solve the assignment problem in a distributed manner. 
It has many variants, and the algorithm originally proposed by Demange, Gale, and Sotomayor [14] (referred to as 
“DGS Auction”) is one of them. Interestingly, DGS Auction algorithm adopts an ascending pricing strategy and 
converges to the minimum price equilibrium [14]. 

In our setting, the transmitters are bidders, and the receivers represent products. Let us look at the assignment 
problem from an auction perspective. The bidders have access to the local channel strength knowledge, i.e., bidder 
i only knows g /C}. Here the left label y^. can be regarded as the price of the product j, meaning that a 

bidder must pay as much as y^. to obtain the product j. For a given price y^., Aij - y^. can be regarded as the 
benefit of the bidder i regarding the project j. We define profit margin by yu^ = maxj{A^j - yvj}^ The objective is 
to deternune the best assignment given this local information, such that each bidder is happy to be assigned to a 
product with the lowest price y^ and in turn highest profit margin yu. Specifically, it is to minimize the price while 
maximizing the profit margin, satisfying yu^ + yvj ^ A^j, Vi, j. 

An algorithm inspired by DGS Auction is detailed in Algorithm 1 where e is a design parameter. If e is small, it 
requires more rounds of iteration to achieve a reasonably “almost optimal” solution. While e is large, as in real 
auctions, the bidder may take risks to pay a non-necessarily high price, leading to a suboptimal solution with a 
faster convergence. A demand set V is maintained among bidders to indicate which bidders are not unassigned any 
product. Oj represents the owner of the product j who successively bids this product. In each round, the bidder 
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who is not assigned any product will bid his most profitable product, i.e., = argmaxj{A^j - yy.} with j* being 

the most profitable product of the bidder i. If the associated profit Aij* - yy.^ is negative or j"" is already assigned 
to the bidder i, then we skip this bidder and consider the next unassigned bidder. Otherwise, bidding process starts. 
If the product was already assigned to another bidder, then this bidder will be added to the demand set and 
reconsidered later. If this product is free, then it will be assigned to the bidder i, and at the same time the price 
of the product will be raised by e. Keep doing this until every bidder has his product without competitors. The 
final assignment of bidders and products is the optimal solution to the assignment problem. 


Algorithm 1 A Decentralized Power Allocation Algorithm via Auction Algorithm 
Require: The bidder i only has the local knowledge A^ = ^an ... an- di ... 
1: Initialization: Set yy. = 0, Oj = 0, Vj, and P = {1, 2 , .. ., A}. 

2: while P ^ 0 do 

3: Choose a bidder i from the demand set, i.e., P ^ P\{i} 

4: For bidder i, find the best values in {Aij - yy. , Vj} 

Wi = max{Aij - yy }, f = arg max{ A^^- - yy } 

3 3 

5: if > e && Oj* ^ i then 

6: if Oj* ^ 0 then 

7: Add i into the demand set, i.e., P ^ P u {Oj*} 

8: end if 

9: Assign the bidder i with the product i.e., Oj* = i 

10: Product raises the price by e, i.e., yy.^ yy.^ + e 

11: end if 

12: end while 


For any set of feasible power allocation parameters [vj = -yuj}j^ there exists an equilibrium price vector 
Yv = {yvj}j- The minimum price equilibrium always exists as long as the corresponding GDoF tuple is feasible. It 
has been proved in [14] that this DCS Auction algorithm leads to the minimum price equilibrium y^, meaning that 
any other feasible price vector is element-wise larger than this minimum equilibrium, i.e., yj^ <yv As such, the 
globally minimal power allocation parameters can be obtained from [r^ = yy.^ - Ai*j*} where (P, j'') belongs to 
the optimal bidder-product assignment. 

If we make e sufficiently small such that the price is raised carefully in each round of auction, [yy. }j obtained 
by Algorithm 1 can achieve arbitrarily close to the minimum equilibrium price. The algorithm converges within a 
finite number of iterations. An adaptive price increasing strategy [38] can be applied to speed up the convergence. 
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TABLE I: System Parameters 


Parameters 

Values 

Cell range 

1km X 1km 

Carrier Frequency 

2.4GHz 

Bandwidth 

5MHz and 10 MHz 

Distance (uniformly distributed) 

[5,30] and [10,60] m 

Transmit Power 

20 and 30 dBm 

Noise power spectral density 

-174 dBm/Hz 

Antenna Height 

1.5 m 

Antenna Gain per Device 

-2.5 dB 

Noise Figure 

7 dB 


VI. Numerical Results 

To demonstrate the gains of our ITLinQ+ mechanism over FlashLinQ and ITLinQ with regard to the sum 
throughput and energy efficiency, we perform numerical analysis under a similar network setup as in [4], [18]. 

As the first setup, we only consider the link scheduling without power control to show the benefit of our new 
decision criterion. We consider a 1 km x 1 km square area and randomly drop n transmitter-receiver pairs. As in 
[4], [18], the simulated channel follows the LoS model in ITU-1411, and the system parameters are listed in Table I. 
Two scenarios are considered: (1) Scenario 1: The distance of any two paired transmitter and receiver is uniformly 
distributed in [5,30]m, all links are supposed to operate over a 5MHz spectrum, and the maximum transmit power 
is 20 dBm; (2) Scenario 2: The range of distance is enlarged to [10,60]m with a larger maximum transmit power 
30 dBm and a wider spectrum bandwidth of lOMHz. We compare ITLinQ+ with no scheduling case where all links 
are activated, FlashLinQ, as well as ITLinQ with properly chosen parameters = 0.7 and M = 25 dB as in [4]. In 
both scenarios, we use the same parameter r] = 0.9 and 7 = 0.1 for our ITLinQ+. 

Fig. 6 shows the sum throughput (bit/sec/Hz) versus the total number of links in the D2D networks ranging from 
8 to 1024. It demonstrates the significant improvement of ITLinQ+ over FlashLinQ and ITLinQ. For instance, with 
in total 1024 links, ITLinQ-r achieves 40% gain in Scenario 1, 60% gain in Scenario 2 over FlashLinQ, and 20% 
gain in Scenario 1 and 40% gain in Scenario 2 over ITLinQ. ^ As a side remark, the design parameters r/, 7 of our 
ITLinQ+ keep unchanged for various scenarios, which meets the demand of robustness in D2D networks. 

To understand better the dramatic energy efficiency improvement of our power control algorithms, we consider a 
smaller random network with 10 link pairs where the desired and cross link strength levels are uniformly distributed 
in [1,2] and [0,1], respectively. We compare in Fig. 7 the sum throughput and energy efficiency (bits per Joule) 
versus SNR. 

Four schemes are compared: (1) Conventional GP-based power allocation (“GP”), (2) distributed power allocation 

^The gap between ITLinQ and ITLinQ+ might be reduced if one can further fine-tune the parameter M in ITLinQ. Unfortunately, there was 
no guideline on how to choose this parameter in [4]. 
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Fig. 6: Sum throughput comparison among no scheduling, FlashLinQ, ITLinQ and our new ITLinQ+ without power control. The 

left is for Scenario 1 and the right for Scenario 2. 




Fig. 7: Sum throughput and energy efficiency versus SNR. 


via DGS-Auction algorithm where the GDoF tuple is found by LP (“LP+AP”), (3) distributed power allocation via 
DGS-Auction algorithm where the GDoF tuple is found by the distributed implementation of GP (“DGP+AP”), and 
(4) full power allocation (“full”) for reference. It is clear that the conventional GP-based power control offers the 
maximal sum throughput, and the assignment-inspired power control with the GDoF tuple obtained by either GP 
or LP gives relatively good performance. Notably, the energy-efficiency of the assignment-inspired power control 
nearly 50 times outperforms that by GP, and 100 times over full power allocation. 

VII. Conclusion 

The GDoF optimality problem of treating interference as noise for Gaussian interference channels has been 
re-formulated in a combinatorial optimization perspective. Thanks to this new formulation, we cast power allocation 
into an assignment problem, which can be solved in polynomial time. A new expression for the TIN-Achievable 
GDoF region is provided, which is more compact and useful than what known before since it eliminates many 
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redundant inequalities. A relaxed version of the condition in [3] on the channel coefficients is given, for which the 
TIN-Achievable GDoF region is a convex poly tope. Finally, a new TIN optimality condition is also revealed, by 
which TIN still achieves the optimal GDoF region for a class of networks different from the one identified in [3]. It 
is also worth noting that our new TIN optimality condition does not violate the conjecture in [3] that the GNAJ 
condition is also necessary “except for a set of channel gain values with measure zero.” 

We are also able to translate these insights into practical communication systems (e.g., D2D networks). With the 
newly found channel strength condition, we employed it as a new decision criterion in a distributed link scheduling 
mechanism. Together with the globally optimal distributed power allocation algorithms, we proposed a distributed 
spectrum sharing mechanism (ITLinQ-r) for D2D networks. By simulation, we have shown that our ITLinQ-i- 
mechanism achieves significant sum throughput improvement over FlashLinQ and ITLinQ with the same level 
implementation complexity. Moreover, ITLinQ-i- also promises a substantial improvement on energy efficiency. 

Appendix 

A. Preliminaries 

1) Weighted Matching: In this work we shall make extensive use of weighted matchings [27] of bipartite graphs. 
We recall here some basic definitions. Let Q = {U^V^£) denote a bipartite graph with left vertices U, right vertices 
V and edges f c ^ x V. A matching TW c is a set of edges, any two of which do not share the same vertex. When 
weights w{u,v) are associated to the edges {u,v) e £, we denote by w(A4) = E(u,v)€M ^(^5 '^) the weight of the 
matching Ai, and we let AT" = arg^ maxu;(Af) denote the maximum weighted matching, i.e., the matching with 
maximum sum-weight. A4* can be characterized as the solution of the integer program: 


max ^ w{u^v)x(u,v), 

{u,v)eS 

(48) 

s.t. ^ x(u, v) < 1, 

ueU:{u,v)€S 

(49) 

E x{u,v)<l, 

veV‘{u,v)eS 

(50) 

x{u,v) e {0,1}. 

(51) 


When equality holds in all constraints, the resulting solution is called a perfect matching, i.e., a matching that covers 
all vertices. The LP relaxation of (48)-(51), obtained by replacing (51) with x{u^v) e [0,1]» is called fractional 
matching [39]. For bipartite graphs, the solution of this LP relaxation is integral, i.e., x g {0,1}. meaning that, given 
a fractional matching, there exists a perfect matching such that the sum-weights of two matchings are equal. In 
other words, there always exists an integral solution to the LP relaxation problem. 

A vertex/edge is called matched if it is involved in a matching; otherwise it is a free vertex/edge. A path is alternating 
if its edges alternate between matched and free edges. The augment operation aug(-) is to exchange matched and 
free edges in an alternating path that starts from and ends to free vertices. For instance, given an alternating path 
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v = {{io,ii),{ii,i 2 ),{i 2 ,h), ■ ■ ■ ^{hnMn+i)) consists of a matching M = {{ii,i 2 ),{h,u), ■ ■ ■ Ai^n-iMn)} and 
free edges V\M, the augment of V results in a new matching M' = aug(P) = 'P\A4 and free edges M. 

2) Geometric Programming: Geometric programing is a powerful tool to solve a class of nonlinear optimization 
problems under a standard form 

min /o(x) 
s.t. /j(x) < 1, i = 

5i(x) = l, z = l,...,p 

where {/i(x), i = 0,1,..., m} are posynomial functions /i(x) : M" h*- R in a form of 

/(xi,... f (53) 

k=l 

with c/c > 0, V/c and {^^(x),7 = 1,... ,p} are monomial functions ^i(x) : M in the form of 

g{xi,...,Xn) = cxl^x^^--x^^ (54) 

with c > 0. 

B. The Kuhn-Munkres Algorithm 

To ease the presentation, we construct a bipartite graph Q = (^,V) with weight of edge (i, j) being and 
ZY, V being transmitter and receiver sets respectively. The Kuhn-Munkres algorithm is to find the maximum weighted 
matching in this bipartite graph. The input is the weight matrix A defined in (17), and the output is the matching 
with maximum sum weights and the corresponding left and right labels {yuj^yvj}j^ in which the left labels {yuj}j 
achieve the maximum left label equilibrium [12], [13]. 

As the initialization, we choose a feasible labeling with yu^ = maxj A^j, V i, and y^. = 0, V j. This labeling 
is feasible, because yu^^ + y^. > Aij always holds for any pair of i eU and j e V. We also construct an equality 
subgraph, Qe, including all the vertices of Q but only those edges (i, j) such that 

yui yvj ~ -^ij 5 Vi G Z7, j € V. (55) 

It has been proved in [12], [13] that, once Qe has a perfect matching M., then this matching Af is a maximum 
weighted matching, and thus the corresponding labels are the final solution to the assignment problem. 

The algorithm consists of multiple rounds of iterations. In each round, we first check if there exists a perfect 
matching in Qe- Note that for a feasible GDoF tuple, the maximum matching always involves edges (j, j) 
according to (11). Thus, the perfect matching consists of edges {(j, j), Vj e 1C}. If not, the left and right labels 
are carefully decreased and increased respectively, and the equality subgraph Qe updated accordingly. Once 
Qe contains the perfect matching {(i,i),V7 ^ the resulting left labels {yuj}j achieve the maximum left 
label equilibrium, yielding the global minimal power allocation Vj = -yuj for all j. The details of the Kuhn- 
Munkres algorithm are given in Algorithm 2, where A/lC^) is the neighborhood of a set of nodes of S in Qe, i.e., 
A/l(5) = {j : (i, j) ^ Qs^^i ^ e V}. An illustrative example is also given as follows. 


(52a) 

(52b) 

(52c) 
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Algorithm 2 A Centralized Power Allocation Algorithm via the Hungarian Method 
Require: Matrix A with ij-th element specified in (17). 

1: Initialization: Set = maxj A^j, V i and y^. = 0, V j. Construct Qe according to {Vui^VvE ^nd choose 

an arbitrary matching M contained in Qe^ Let S = T = 0. 

2: \i3 M = Vj e K,} in Qe then 

3: Ti = -yuii V i, and return 

4: else 

5: Pick a free vertex ueU 

6: S ^ 

7: end if 

8: if = T then = mmi^sjtriVui + Vvj - Mj), 

9 : Update Uu^ ^ yu^ - ai, if keS, 

10 : Update y^^ + at, if k^T, 

11: Update M = {(t,j) : Vm + Vv^ = Ay} 

12: else 

13: Pick V G AfL{S)\T 

14: if is a free vertex then 

15: Augment the alternating path v that contains the matching M 

16: update Ai ^ aug({it ^ 'i^}) goto 2 

17: else 

18: if is matched to u' then 

19: 5 ^ 5 u {u'}, T {'i;} 

20: goto 8 

21 : end if 

22: end if 

23: end if 


Example 3. The detailed power allocation procedure according to Algorithm 2 is presented as follows. 

As an initialization, we assign 

Vui = 1 - 5 , yu 2 = 0 - 5 , Vus = Vv^ = yV 2 =yV 3 =^ ( 56 ) 

such that we construct the equality subgraph Qe with edges {(1,1), (2,3), (3,1)} as in Fig. 8(b). Note that Qe 
does not contain a perfect matching. So we choose an arbitrary matching, e.g., M. = {(2,3), (3,1)}, as shown in 
Fig. 8(b). 

In the first round, we choose a free vertex ui in Qe and set S = {i^i} and T = 0. Because A/l( 5) = {i;i} ^ T, we 
go to line-12 and pick vi e A/l( 5)\T. Note that vi is matched to us, and thus we update S = {ui^us}, T = 
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Fig. 8: Illustration of the Hungarian method, (a) The weighted bipartite graph with weight of zj-th edge being Aij. (b) The 
equality subgraph with initiated labeling, (c) The updated equality subgraph with first round label update, (d) The update 
equality subgraph with second round label update. The edges in red give a matching in equality subgraphs. 


As A/l( 5) = = T, we go to line-8 and obtain 

Q!L = min - Ajj} = 0.2 (57) 

Uie{ui,U3},Vje{v2,V3} 

As such, we have 

Vui = 1-3, yu2 = 0.5, Vus = 0.8, = 0.2, = 0, y^^ = 0. (58) 

and go to line-2. 

In the second round, we update the equality subgraph Qe with edges {(1,1),(2,3),(3,1),(3,3)} and there still 
does not contain a perfect matching, as shown in Fig. 8(c). Thus, we choose M. = {(2,3), (3,1)} as a matching. Still, 
we pick the free vertex ui, and set S = {txi} and T = 0. Now, as ^ T, we pick vi e A/l(‘5)\T. 

Because vi is matched to its, we update S = {ui,us}, T = At this point, A/l(5) = {vi^vs} ^ T again. 
Thus, we pick vs e A/l(‘ 5)\T. Due to vs is matched to U 2 , we update S = {iti,it 2 ,its} and T = {i’i,t;s}. Here 
A/l( 5) = {t;i,i’s} = T, we go to line-8 and have 

aL= ^ min - Ai^} = 0.1 (59) 

Ui€{ui ,U2 ,U3} ,Vj e{v2} 

As such, the labels are updated as 

Vui = 1-2, = 0.4, 2/„3 = 0.7, = 0.3, = 0, y„^ = 0.1. (60) 

and then we go to line-2. 

Till now, in the updated equality subgraph shown in Fig. 8(d), we have a perfect matching {(1,1), (2,2), (3,3)}. 
Thus, the algorithm returns with 

n = -1.2, r2 = -0.4, rs = -0.7. (61) 


C. Proof of Theorem 3 

For the sake of this proof, we denote by Vs the region defined by (5), and by vj^^^ the region defined by (21). 
Our goal is to show that Vs = Vj^^^ for any S 
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Vs ^ To prove this, we show that for any inequality presented in we can always find the same 

one in Vs^ 

Given a subnetwork ^[5], the matching with the maximum weight is a degree-1 subgraph. ^ Connecting the 
direct links will lead to single or multiple disjoint cycles, as all the nodes has degree-2. 

For the single-cycle case, this cycle corresponds to a sum GDoF constraint in Vs^ For the multiple-cycle case, 
each cycle corresponds to a sum GDoF constraint in Vs of the users involved in this cycle. Thus, the sum GDoF 
constraint with the maximum weighted matching in corresponds to the combination of these sum GDoF 

constraints in Vs^ 

As such. Vs contains or implies all the constraints in As Vs has more constraints, it follows that 

Vs ^ Vf^^. 

pTiNA g prove this, we show that, for any subset of users 5, the TINA GDoF region confined by Vs is 

no larger than that by it is clear that Vs is determined by individual GDoF and sum GDoF constraints of 

any subset of users in S. The individual GDoF constraints of two regions are identical. Thus, our focus will be on 
the sum GDoF constraints for users in S with |5| > 2. 

For the user set 5, the sum GDoF constraints in Vs only come from (1) the sum GDoF constraints with all 
possible permutations of 5, and (2) the combination of a number of individual and/or sum GDoF constraints of 
subsets of S. For the first case, the sum GDoF constraint in Vs is dominated by the maximum weight of any possible 
matchings (associated with cyclic sequences). For the second case, suppose the combination involves a number of 
subnetworks, where the subnetworks may have intersections. This combination of constraints involves every user 
with equal times (say b times), otherwise, the combination will not lead to a sum GDoF constraint, because it is a 
weighted sum GDoF constraint and can be implied by the combination of other sum GDoF constraints. Each sum 
GDoF constraint for a subnetwork involves a cyclic sequence and hence forms a matching. Thus, the combination 
of sum GDoF constraints corresponds to a fractional perfect matching by assigning x{u,v) in (48)-(51) with In 
bipartite graphs, the weight of any fractional perfect matching equals the weight of a perfect matching [27], [39]. 

Thus, neither the weight of any matching nor of any fractional matching is greater than the maximum weighted 
matching, such that the sum GDoF constraints in will be more restrictive than those or any combinations in 

Vs, i.e., V^^^^ c This completes the proof. 

D. Proof of Theorem 4 

In what follows, we prove that under condition (23), is monotonically increasing. Hence, from (6) this 

immediately implies that inspection, is a convex poly tope. 

Let us start with \S\ = 2. Suppose without loss of generality S = {k,j}. Due to the condition (23), mm{akk,<^jj} > 
akj + CYjk, then it is easy to verify that and Vj^^^ ^ 

Then, we prove the general cases with the following lemma. 

^If direct links are in the matching, we can eliminate them from S, which does not affect our proof. 




31 


Lemma 1. Given a subgraph ^[5] with weights g 5}, the difference of maximum weighted matching with 

and without the user k is bounded by 

w{M*s) - < . max {an, + auj - a[j} 

Proof: Suppose without loss of generality that the maximum weighted matching of ^[5] {k e S) includes links 
(i, /c) and {k^j) with weights and a'^- respectively and i^jtk. Note that whether i = j or not does not affect 
our proof. After removing user k and edges {k,j) from the matching, and adding the link (i, j) with weight 

we have a matching for S\{k}. Thus, for all {(i,k),{k,j)} e we have 




< max {a^k + Q(L- - oi'jA 

(62) 

— . \^^ik ^kj ~ ^ij } 

(63) 

= max {ttife + akj - a' ■}. 

i,j£S ,i,j^k 

(64) 


■ 

Together with the condition (23), we have 


akk > w{M*s) - w{M*s\{k}) 

(65) 


for any user k eS. Thus, by comparing the sum GDoF constraints without and with user k 


j€S\{k} j^S\{k} 

E ^ E “ii - (67) 

j^<s j€S 

= E ajj +akk-w{M*s), (68) 

jeS\{k} 

it is readily verified that as long as (23) is satisfied, the sum GDoF constraint with user k is not implied by the 
sum GDoF constraint without user k and the individual GDoF constraint dk < oikk- In other words, with user k, the 
GDoF region is not decreasing. It follows immediately that ^ {\fk e S). More generally, if Si Q S 2 , 

then c This completes the proof. 

E. Proof of Theorem 5 

Due to the fact that IV 2 and that, under condition (23), = 7^;^^^^, achievability trivially follows. 

For the converse, we follow the cyclic outer bounds first revealed in [40, Theorem 2] and later used to prove the 
optimality of TIN condition in [3, Theorem 3]. 
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Thus, for the if-user Gaussian interference channel in the weak interference regime, the GDoF region under the 
condition (23) is included in the set of GDoF tuples (di, (i 2 ,..., such that 


dj<ajj, Vjf G/C (69) 

E ^ min{/7r 5 ^7r,0: • • • 5 9n,m-l }, (70) 

i=o 

for any ordered subset tt = (io,H, • • • where we define 

m-l 

/^ = E (71) 

i=o 

m-l 

97T,k = E ~ > fc = 0,...,TO-l, (72) 

i=0 


and where the index subscript arithmetic is modulo m. 

When m = 2, then condition (23) is equivalent to the GNAJ condition and the bound is known to be tight. When 
m > 2, let us first consider the bound formed by the “g” terms in (70). Notice that the left-hand side of (70) depends 
only on the indices in tt but not on its order. Hence, letting S denote a given unordered subset of size m of /C 
and using the short-cut notation tt g 7r(5) to indicate the ordered sets formed with the elements of 5, i.e., the 
permutations of 5, we can write 


min min { Qtt a- I 

7T€7T{S)k=0,...,m-l^ ’ ^ 


= E 


a 


JJ 


[ m-l 



E 

[j=o 

j ^ij-i 


1 


■lij-, min 

- max \ 

E 

7re7r(S) | 

[i=o 





(73) 


where (73) is due to the condition (24). If the maximum weighted matching involves multiple cycles, then the sum 
GDoF outer bound can be the combination of multiple sum GDoF constraints associated with the corresponding 
cyclic sequences. Thus, (73) still holds, because condition (24) holds for any subset of S ^ 1C. Due to the fact that 

m-l m-l 

S ~ Clijij+i 5 ^ijij ~ } 

j=0 j=0 


we have that 


Y,dj< min min{/7,,^7 ,,o,-..,^7 r,m-i} (74) 

j^S 7ren{S) 


= (75) 

which coincides with for every S ^ 1C. Under the condition (23), the TINA is the largest polyhedral region, 

so the converse bound is tight. 
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